Identification of a 3-Gene Prognostic Index for Papillary Thyroid Carcinoma.

The accurate determination of the risk of cancer recurrence is a critical unmet need in managing thyroid cancer (TC). Although numerous studies have successfully demonstrated the use of high throughput molecular diagnostics in TC prediction, it has not been successfully applied in routine clinical use, particularly in Chinese patients. In our study, we objective to screen for characteristic genes specific to PTC and establish an accurate model for diagnosis and prognostic evaluation of PTC. We screen the differentially expressed genes by Python 3.6 in The Cancer Genome Atlas (TCGA) database. We discovered a three-gene signature Gap junction protein beta 4 (GJB4), Ripply transcriptional repressor 3 (RIPPLY3), and Adrenoceptor alpha 1B (ADRA1B) that had a statistically significant difference. Then we used Gene Expression Omnibus (GEO) database to establish a diagnostic and prognostic model to verify the three-gene signature. For experimental validation, immunohistochemistry in tissue microarrays showed that thyroid samples' proteins expressed by this three-gene are differentially expressed. Our protocol discovered a robust three-gene signature that can distinguish prognosis, which will have daily clinical application.

Network Pharmacology to Explore the Molecular Mechanisms of Prunella vulgaris for Treating Hashimoto's Thyroiditis.

Purpose: Prunella vulgaris (PV), a traditional Chinese medicine, has been used to treat patients with thyroid disease for centuries in China. The purpose of the present study was to investigate its bioactive ingredients and mechanisms against Hashimoto's thyroiditis (HT) using network pharmacology and molecular docking technology to provide some basis for experimental research. Methods: Ingredients of the PV formula were retrieved from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. Additionally, HT-related genes were retrieved from the UniProt and GeneCards databases. Cytoscape constructed networks for visualization. A protein-protein interaction (PPI) network analysis was constructed, and a PPI network was built using the Search Tool for the Retrieval of Interacting Genes (STRING) database. These key targets of PV were enriched and analyzed by molecular docking verification, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Results: The compound-target network included 11 compounds and 66 target genes. Key targets contained Jun proto-oncogene (JUN), hsp90aa1.1 (AKI), mitogen-activated protein kinase 1 (MAPK1), and tumor protein p53 (TP53). The main pathways included the AGE-RAGE signaling pathway, the TNF signaling pathway, the PI3K-Akt signaling pathway, and the mitogen-activated protein kinase signaling pathway. The molecular docking results revealed that the main compound identified in the Prunella vulgaris was luteolin, followed by kaempferol, which had a strong affinity for HT. Conclusion: Molecular docking studies indicated that luteolin and kaempferol were bioactive compounds of PV and might play an essential role in treating HT by regulating multiple signaling pathways.

Prioritizing Susceptible Genes for Thyroid Cancer Based on Gene Interaction Network.

Thyroid cancer ranks second in the incidence rate of endocrine malignant cancer. Thyroid cancer is usually asymptomatic at the initial stage, which makes patients easily miss the early treatment time. Combining genetic testing with imaging can greatly improve the diagnostic efficiency of thyroid cancer. Researchers have discovered many genes related to thyroid cancer. However, the effects of these genes on thyroid cancer are different. We hypothesize that there is a stronger interaction between the core genes that cause thyroid cancer. Based on this hypothesis, we constructed an interaction network of thyroid cancer-related genes. We traversed the network through random walks, and sorted thyroid cancer-related genes through ADNN which is fusion of Adaboost and deep neural network (DNN). In addition, we discovered more thyroid cancer-related genes by ADNN. In order to verify the accuracy of ADNN, we conducted a fivefold cross-validation. ADNN achieved AUC of 0.85 and AUPR of 0.81, which are more accurate than other methods.

Significance of DMBT1 in Papillary Thyroid Carcinoma Concurrent With Hashimoto's Thyroiditis.

BACKGROUND: Papillary thyroid carcinoma (PTC) concurrent with Hashimoto's thyroiditis (HT) was associated with a better clinical prognosis. This study aimed to investigate a potential mRNA gene that affects the development of PTC, which helps PTC concurrent with HT patients have a better prognosis. MATERIAL/METHODS: PTC data were obtained from The Cancer Genome Atlas (TCGA) database. And the validation data of tissue specimens were collected from Guangzhou First People's Hospital. The thyroid tissue sections were hybridized with deleted in malignant brain tumor 1 (DMBT1) probes by situ hybridization. Survival rates were analyzed using Kaplan-Meier curves, and the log-rank test was used to compare group survival rates. Prognosis clinicopathological factors were analyzed by Cox regression. Gene Ontology (GO) and Kyoto Gene and Genomic Encyclopedia (KEGG) pathway enrichment analyses were performed using single-sample gene set enrichment analysis (ssGSEA). Finally, the correlation of deletion in DMBT1 expression with overall immune status, tumor purity, and human leukocyte antigen (HLA) gene expression profile was analyzed. RESULTS: HT was significantly associated with sex, tumor foci, extrathyroidal extension (ETE), residual tumor, and tumor stage (T stage). Moreover, PTC concurrent with HT had a lower risk of recurrence versus non-HT groups. A total of 136 differentially expressed mRNAs (DEMs) were identified between HT and non-HT groups. Among them, the expression level of DMBT1 in HT groups was statistically higher than that in non-HT groups. A significant association with ETE and recurrence was revealed in the high expression and the low expression of DMBT1. Furthermore, DMBT1 was an independent predictor of survival. The overall immune activity of high expression of DMBT1 was higher than that of the low-expression group. CONCLUSIONS: The PTC patients with HT had better behavior features and prognosis than those with simple PTC. DMBT1 in PTC-HT patients was a potential possible factor that inhibits tumors. High expression of DMBT1 may improve PTC prognosis by immune-related pathways.

Prediction of Genetic Factors of Hyperthyroidism Based on Gene Interaction Network.

The number of hyperthyroidism patients is increasing these years. As a disease that can lead to cardiovascular disease, it brings great potential health risks to humans. Since hyperthyroidism can induce the occurrence of many diseases, studying its genetic factors will promote the early diagnosis and treatment of hyperthyroidism and its related diseases. Previous studies have used genome-wide association analysis (GWAS) to identify genes related to hyperthyroidism. However, these studies only identify significant sites related to the disease from a statistical point of view and ignore the complex regulation relationship between genes. In addition, mutation is not the only genetic factor of causing hyperthyroidism. Identifying hyperthyroidism-related genes from gene interactions would help researchers discover the disease mechanism. In this paper, we purposed a novel machine learning method for identifying hyperthyroidism-related genes based on gene interaction network. The method, which is called "RW-RVM," is a combination of Random Walk (RW) and Relevance Vector Machines (RVM). RW was implemented to encode the gene interaction network. The features of genes were the regulation relationship between genes and non-coding RNAs. Finally, multiple RVMs were applied to identify hyperthyroidism-related genes. The result of 10-cross validation shows that the area under the receiver operating characteristic curve (AUC) of our method reached 0.9, and area under the precision-recall curve (AUPR) was 0.87. Seventy-eight novel genes were found to be related to hyperthyroidism. We investigated two genes of these novel genes with existing literature, which proved the accuracy of our result and method.

Tumor Mutation Burden Predicts Relapse in Papillary Thyroid Carcinoma With Changes in Genes and Immune Microenvironment.

Background: The risk factors of papillary thyroid carcinoma (PTC) recurrence are meaningful for patients and clinicians. Tumor mutation burden (TMB) has been a biomarker for the effectiveness of immune checkpoint inhibitor (ICI) and prognosis in cancer. However, the role of TMB and its latent significance with immune cell infiltration in PTC are still unclear. Herein, we aimed to explore the effect of TMB on PTC prognosis. Material and Methods: RNA-seq and DNA-seq datasets of PTC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The Gene Ontology (GO) and gene set enrichment analysis (GSEA 4.0.1) were applied further to explore potential differences in PTC patients' biological functions. The differentially expressed genes (DEGs) and immune microenvironment between the high and low TMB groups were determined. Results: TMB had the highest AUC score than other clinical indicators in ROC analysis on recurrence-free survival, and a higher TMB score was related to a worse prognosis. Further, GSEA showed a higher level of oxidative phosphorylation (OXPHOS) in the high TMB group, and four genes correlated with recurrence-free survival rate were identified. The abundance of CD8+ T cells and M1 macrophages in the high TMB group was significantly lower than that in the low TMB group. Conclusions: Our study found that TMB was a better predictor variable at evaluating the risk of PTC recurrence. Moreover, TMB-related genes conferred dramatically correlated prognosis, which was worth exploring in guiding postoperative follow-up and predicting recurrence for PTC patients.

Identifying Thyroid Carcinoma-Related Genes by Integrating GWAS and eQTL Data.

Thyroid carcinoma (TC) is the most common endocrine malignancy. The incidence rate of thyroid cancer has increased rapidly in recent years. The occurrence and development of thyroid cancers are highly related to the massive genetic and epigenetic changes. Therefore, it is essential to explore the mechanism of thyroid cancer pathogenesis. Genome-Wide Association Studies (GWAS) have been widely used in various diseases. Researchers have found multiple single nucleotide polymorphisms (SNPs) are significantly related to TC. However, the biological mechanism of these SNPs is still unknown. In this paper, we used one GWAS dataset and two eQTL datasets, and integrated GWAS with expression quantitative trait loci (eQTL) in both thyroid and blood to explore the mechanism of mutations and causal genes of thyroid cancer. Finally, we found rs1912998 regulates the expression of IGFALS (P = 1.70E-06) and HAGH (P = 5.08E-07) in thyroid, which is significantly related to thyroid cancer. In addition, KEGG shows that these genes participate in multiple thyroid cancer-related pathways.

ALDH5A1 acts as a tumour promoter and has a prognostic impact in papillary thyroid carcinoma.

Thyroid cancer is the most common endocrine carcinoma, with papillary thyroid carcinoma (PTC) accounting for 80%-90% of thyroid cancers. Accumulating studies reported that mitochondria plays an important role in the regulation of cell proliferation. ALDH5A1, may function as an oncogene or tumour suppressor in various human cancers, and the role of ALDH5A1 in PTC is still unclear. The aim of this study was to investigate the clinical significance of ALDH5A1 expression and its functions in PTC. In this present study, we studied ALDH5A1 expression on primary papillary thyroid carcinoma (PTC) in The Cancer Genome Atlas (TCGA) database. Results showed that the levels of ALDH5A1 were found positively correlated with tumour stage, metastasis, lymph node stage, and higher levels of ALDH5A1 demonstrated poor disease-free survival (DFS). Immunohistochemistry (IHC) revealed that significantly higher expression of ALDH5A1 was found in PTC tissues. On the other hand, knockdown of ALDH5A1 significantly inhibited PTC cell proliferation, migration and invasion detection found the migration and invasion of cells also were hindered when ALDH5A1 level was reduced. The knockdown of ALDH5A1 inhibited the expression of Vimentin and promoted the expression of E-cadherin. In brief, knockdown of ALDH5A1may act as a novel molecular target for the prevention and treatment of PTC. SIGNIFICANCE OF THE STUDY: The present study focused on the role and the potential mechanism of ALDH5A1 in papillary thyroid carcinoma. We demonstrated that reduced expression of ALDH5A1 might inhibit the progression of TC by inhibiting cell proliferation, migration and invasion and reversing epithelial-mesenchymal transition (EMT). The findings ensured the interaction relation between ALDH5A1 and EMT in PTC, providing a novel biological marker for PTC and enriching the potential strategies for TC treatment.

MicroRNA-3690 promotes cell proliferation and cell cycle progression by altering DKK3 expression in human thyroid cancer.

An increasing amount of evidence has demonstrated the importance of microRNAs (miRNAs/miRs) in the tumorigenesis of malignant types of cancer, and data retrieved from The Cancer Genome Atlas database revealed that miR-3690 was upregulated in thyroid cancer (TC). The present study focused on the biological function and mechanism of miR-3690 in TC, demonstrating that miR-3690 expression was significantly elevated in TC cells and clinical tissues. Functional studies indicated that miR-3690 acted as an oncogene in TC by promoting cell proliferation, colony formation and cell cycle progression in association with the increased expression of cyclin E and c-myc. Mechanistically, prediction software indicated that Dickkopf-related protein 3 (DKK3) was a target of miR-3690, which was confirmed by the results of luciferase reporter assays and western blotting. DKK3 silencing abrogated the functions of miR-3690-in on TC cell proliferation. Collectively, the findings of the present study demonstrated that miR-3690 promoted TC cell proliferation and indicated miR-3690 as a potential biomarker and therapeutic target for TC.

Potential five-mRNA signature model for the prediction of prognosis in patients with papillary thyroid carcinoma.

Although the mortality rate of papillary thyroid carcinoma (PTC) is relatively low, the recurrence rates of PTC remain high. The high recurrence rates are related to the difficulties in treatment. Gene expression profiles has provided novel insights into potential therapeutic targets and molecular biomarkers of PTC. The aim of the present study was to identify mRNA signatures which may categorize PTCs into high-and low-risk subgroups and aid with the predictions for prognoses. The mRNA expression profiles of PTC and normal thyroid tissue samples were obtained from The Cancer Genome Atlas (TCGA) database. Differentially expressed mRNAs were identified using the 'EdgeR' software package. Gene signatures associated with the overall survival of PTC were selected, and enrichment analysis was performed to explore the biological pathways and functions of the prognostic mRNAs using the Database for Visualization, Annotation and Integration Discovery. A signature model was established to investigate a specific and robust risk stratification for PTC. A total of 1,085 differentially expressed mRNAs were identified between the PTC and normal thyroid tissue samples. Among them, 361 mRNAs were associated with overall survival (P<0.05). A 5-mRNA prognostic signature for PTC (ADRA1B, RIPPLY3, PCOLCE, TEKT1 and SALL3) was identified to classify the patients into high-and low-risk subgroups. These prognostic mRNAs were enriched in Gene Ontology terms such as 'calcium ion binding', 'enzyme inhibitor activity', 'carbohydrate binding', 'transcriptional activator activity', 'RNA polymerase II core promoter proximal region sequence-specific binding' and 'glutathione transferase activity', and Kyoto Encyclopedia of Genes and Genomes signaling pathways such as 'pertussis', 'ascorbate and aldarate metabolism', 'systemic lupus erythematosus', 'drug metabolism-cytochrome P450 and 'complement and coagulation cascades'. The 5-mRNA signature model may be useful during consultations with patients with PTC to improve the prediction of their prognosis. In addition, the prognostic signature identified in the present study may reveal novel therapeutic targets for patients with PTC.

Targeting of TLE3 by miR-3677 in human breast cancer promotes cell proliferation, migration and invasion.

Numerous studies have indicated an important function of microRNAs (miRs) in breast cancer (BC) progression, oncogenesis and metastasis. However, the function of miR-3677, which has been revealed to be upregulated in BC [The Cancer Genome Atlas (TCGA) data], has not been investigated to date. In the present study, miR-3677 was revealed to be upregulated in BC as determined using TCGA. miR-3677 was significantly upregulated in BC tissues and cell lines compared with those noted in adjacent non-cancerous tissues and primary normal breast cells (P<0.05). The overexpression of miR-3677 promoted the cell proliferation, migration and invasion of BC cells. Using bioinformatics algorithms and luciferase assays, a novel target gene for miR-3677, namely transducin-like enhancer of Split3 (TLE3), was identified. Silencing of TLE3 in miR-3677-transfected BC cells suppressed their proliferation and migration. An inverse correlation was observed between miR-3677 and TLE3 expression levels in human BC tissues. In conclusion, the present study demonstrated that miR-3677 promoted BC cell proliferation, migration and invasion by inhibiting TLE3 expression, which provided a novel mechanism and a promising therapeutic target for patients with BC.

Prognostic implications of the BRAF-V600E mutation in papillary thyroid carcinoma based on a new cut-off age stratification.

The BRAF-V600E mutation is the most common and specific oncogenic event known in papillary thyroid carcinoma (PTC). However, it remains controversial whether there is an association between the BRAF-V600E mutation and clinicopathologically aggressive characteristics of PTC. The purpose of the present retrospective study was to investigate the significance of the BRAF-V600E mutation in predicting prognostic and aggressive clinicopathological characteristics according to a new age-based stratification. Clinical data and the BRAF-V600E mutation status of 475 patients with PTC were downloaded from The Cancer Genome Atlas database. The association between BRAF-V600E status and clinicopathological characteristics was analyzed by χ2 test or Fisher's exact test. Recurrence-free survival rate (RFS) was analyzed using the Kaplan-Meier method. Aggressive clinicopathological factors associated with recurrence were analyzed by Cox multivariate regression. This study was conducted on 219 cases of patients with PTC with a known BRAF-V600E mutational status. In the ≥55 years age group, BRAF-V600E was found to be significantly associated with aggressive PTC characteristics, including tumor size, PTC subtype, radioactive iodine (RAI) dose, follow-up time, recurrence, recurrence risk stage, advanced T stage, advanced N stage and American Joint Committee on Cancer (III/IV) stage (all P<0.05). RFS was analyzed by the log-rank test and exhibited statistically significant differences in the ≥55 years group (P=0.041), but there was no significant difference in the <55 group (P=0.757), according to the BRAF-V600E mutation status. The BRAF-V600E gene was excluded from the recurrence Cox multivariate regression model. The BRAF-V600E mutation was found to better predict aggressive and recurrent PTC based on age stratification with the cut-off age of 55 years. The synergistic interaction between BRAF-V600E mutation and the new age stratification may help clinicians reach the optimal decision in terms of surgical approach and the extent of surgery.

Survival of aggressive variants of papillary thyroid carcinoma in patients under 55 years old: a SEER population-based retrospective analysis.

BACKGROUND: Patients younger than 55 years of age with papillary thyroid carcinoma (PTC) have excellent survival. Diffuse sclerosing variant (DSV) and tall cell variant (TCV) of PTC are associated with aggressiveness; the survival of patients <55 years of age with these variants is still unclear. We aim to investigate the clinicopathological features and survival of these variants in the age group <55 years. METHODS: All adult patients (<55 years old) with DSV, TCV and conventional PTC (CPTC) came from the Surveillance, Epidemiology, and End Results program (1988-2013). Kaplan-Meier method and log-rank test were used to analyze the survival. Prognostic factors associated with survival were analyzed by Cox multivariate regression. RESULTS: There were 280 DSV, 615 TCV, and 56287 CPTC in the age group <55 years. DSV and TCV were associated with multifocality, extrathyroidal extension, lymph node and distant metastasis (all p < 0.05). The 10-year disease-specific survival (DSS) of TCV was worse than CPTC (96.3 vs. 99.4%, p < 0.01), but there was no significant difference between DSV and CPTC (99.5 vs. 99.4%, p > 0.05). Cox multivariate regression showed TCV was the independent predictor of DSS (HR: 5.39, p < 0.01). CONCLUSION: In the age group <55 years, DSV and TCV are more likely to exhibit aggressive characteristics than CPTC. Patient <55 years of age with DSV have excellent survival likewise, while patients <55 years of age with TCV carry worse survival. Further investigation for the recurrence risk of patients <55 years with these variants would contribute to optimal clinical management making.

The role of two tumor foci for predicting central lymph node metastasis in papillary thyroid carcinoma: A meta-analysis.

BACKGROUND AND OBJECTIVE: Two tumor foci are the most common in multifocal papillary thyroid carcinoma, but whether they should be regarded as the indicator of central lymph node metastasis remains unclear. To investigate the role of two tumor foci for predicting central lymph node metastasis, we performed a meta-analysis of published studies. METHODS: We performed a systematic literature search of PubMed, Embase and Web of Science prior to September 29, 2017. The relevant articles were examined and the eligible studies were included to assess the metastatic risk of central lymph node in papillary thyroid carcinoma with one, two and more than two (>2) tumor foci. RESULTS: Five eligible studies included 4045 patients in this meta-analysis. Two tumor foci were the most common in multifocal papillary thyroid carcinoma (63.8%, 939/1471). Multifocality group showed a higher risk of central lymph node metastasis compared with unifocality group (odds ratio: 1.58, 95% confidence interval: 1.37-1.81). The risk of central lymph node metastasis was higher in two tumor foci group than unifocality group (odds ratio: 1.38, 95% confidence interval: 1.17-1.62). However, this risk in two tumor foci group was lower than >2 tumor foci group (odds ratio: 0.62, 95% confidence interval: 0.42-0.92). Begg's test revealed no obvious publication bias. CONCLUSION: This meta-analysis suggested that two tumor foci should be regarded as the predictive factor of central lymph node metastasis, but the role of it was less important than three or more than three tumor foci. Understanding the role of two tumor foci for predicting central lymph node metastasis may help clinicians make an optimal decision of treatment and the extent of surgery for multifocal papillary thyroid carcinoma.